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RESUMEN 

Presentamos las eficiencias de pcrdida de masa y mctales obtenidas de mode- 
los numericos de vientos galacticos empujados por la energia depositada en brotcs de 
formacion estelar nucleares. Los brotes de formacion estelar contienen ciimulos es- 
telares jovenes los cuales inyectan la energia suficiente para empujar parte del medio 
interestelar fuera de las galaxias. En algunos casos los vientos galacticos contienen 
una importante parte de los metales producidos por las nuevas generaciones este- 
lares. Para estudiar las eficiencias de perdida de masa y metales hemos desarrollado 
simulaciones numercias 3D "N-Cuerpos/Smooth Particle Hydrodynamics" de vien- 
tos galacticos para los casos: adiabaticos y con perdidas radiativas. Los modelos 
numericos cubren una amplio intervalo de masas en los brotes de formacion estelar 
(de ~10^ a ~10^ M©) y de masas en las galaxias anfitrionas (de ^6x10^ a ^10^^ 
M©). Las regiones de formacion estelar concentradas en el centro del potential, son 
una maquinaria importante para la pcrdida y redistribucion de masa y metales en 
este tipo de galaxias. 

ABSTRACT 

We present results from models of galactic winds driven by energy injected by 
nuclear starbursts. The total energy of the starburst is provided by young central 
stellar clusters and parts of the galactic interstellar medium are pushed out as 
part of the galactic wind (in some cases the galactic wind contains an important 
part of the metals produced in the new generation of stars). We have performed 
adiabatic and radiative 3D N-Body /Smooth Particle Hydrodynamics simulations 
of galactic winds using the GADGET-2 code. The numerical models cover a wide 
range of starburst (from ^10^ to ~^10^ M©) and galactic gas masses (from ^--^6x10^ 
to ~10^^ Mq). The concentrated central starburst regions are an efficient engine 
for producing of the mass and metal loss in galaxies, and also for driving the metal 
redistribution in the galaxies. 
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1. INTRODUCTION 

A great amount of evidence both, observa- 
tional and theoretical, indicates that Galactic Winds 
(GWs) are a necessary and very important ingredi- 
ent of the evolution of galaxies and the inter-galactic 
medium. Their presence can explain, among other 
things, the very low metallicities found in dwarf 
galaxies or the abundance of metals observed in the 
intra-cluster medium. Such outflows have been ob- 



'^Instituto de Ciencias Nucleares, Universidad Nacional 
Autonoma de Mexico 

^Centro de Radioastronomia y Astrofisica, Universidad 
Nacional Autonoma de Mexico. 



served in galaxies of a variety of masses and at differ- 
ent redshifts (e.g., Martin 1999, Pettini et al. 2001, 
Veilleux et al. 2005). At the same time, different 
analytical and numerical models have been put for- 
ward to explain these outflows (e.g.. Chevalier & 
Clegg 1985, Tomisaka & Ikeuchi 1988, Heckman et 
al. 1990, De Young & Heckman 1994, MacLow & 
Ferrara 1999, Strickland & Stevens 2000, Canto et 
al. 2000, Tenorio-Tagle et al. 2003, Fujita 2004, 
Cooper et al. 2008, Fujita et al. 2009, Rodriguez- 
Gonzalez et al. 2009). 

Long ago, Lynds & Sandage (1963) found a 
large outflow in M82 and gave the first definition 
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of a starburst (SB). They proposed that material 
was ejected from the nuclear regions as a conse- 
quence of a series of SN explosions. Since then, 
galaxies with concentrated central star bursts (M82, 
NGC 253, NGC 1569, to name a few) have been used 
to study the GW phenomenon. 

Tcnorio-Taglc ct al. (1998) showed the impor- 
tance of the galactic potential component due to the 
cold dark matter halo in realistic models of GW, and 
concluded that the fate of the stellar wind and SNs 
material depends on the galactic properties and on 
the total energy injected by the SB. They showed 
that dwarf galaxies with an interstellar medium 
(ISM) with masses of the order of 10^ M©, retain 
their metals unless they undergo an extreme burst 
of star formation, much larger than those presently 
observed. On the other hand, MacLow & Ferrara 
(1999) explored both, analytically and numerically, 
the effects of stellar winds and SN explosions, in SBs 
with mechanical luminosities in the '--^10'^^ - 10'^^ 
erg s~^ range, on the ISM of dwarf galaxies with 
gas masses in the ^ 10^ to 10^ Mq range. They 
distinguished between two possibilities of gas ejec- 
tion: a) A "blowout" , which require material be ac- 
celerated above the escape speed at a distance of 
about three times the scale height, b) A "blow- 
away" , in which virtually all the ISM is accelerated 
and lost to the galaxy. The blowout would hap- 
pen preferably for flatter galaxies, rounder galaxies 
are more likely to be completely disrupted (blow- 
away) if enough energy is injected. Their conclusion 
was that dwarf galaxies in that mass range undergo 
blowout for modcratc-to-high luminosity starbursts 
values {L ^ 10'^®- 10'^^ erg s~^), whereas for galax- 
ies with gas masses Alg > 10^ Mq and mechanical 
luminosity L ^ 10"^^ erg s~^ the blowout was in- 
hibited. They also found that in the lowest mass 
objects, blow-away occurs virtually independently of 
L. 

Many authors have shown that in the nearest 
starburst galaxies GWs could form through the col- 
lective effect of many individual stellar cluster winds, 
which in turn are formed by the collective effect of 
many individual stellar winds. Perhaps, the best 
studied starburst galaxy is M82, which has an ex- 
tended, biconical filamentary structure in the optical 
(Shopbell & Bland-Hawthorn 1998; Ohyama et al. 
2002). In M82, the super stellar clusters (SSCs) Ha 
luminosities are in the (0.01-23) x 10"^^ erg s~^ range 
(Melo et al 2005). These authors compiled a cata- 
log of 197 SSCs in M82 and 48 SSCs in NGC 253 
(Melo 2005), with masses in the lO^* < M/M^ < 10^ 
range (more recently, additional young star clusters 



in M82 have been reported by Smith et al. 2007 
and Westmoquette et al. 2007). Other examples 
are NGC 1569 (Anders et al. 2004; Westmoquette, 
Smith & Gallagher HI 2008) and M83 (Harris et 
al. 2001), with SSCs densities and masses similar to 
those of M82. Such an exceptional density of mas- 
sive clusters (i.e '--^ 620 kpc~^, see also Anders et al 
2004; Melo et al. 2005; Melo 2005) is the best sce- 
nario to study powerful GWs driven by SSC winds 
originated in nuclear starburst regions. 

Yet, another interesting effect is when the mate- 
rial ejected from the SSCs winds does not abandon 
the galaxy. Tenorio-Tagle 1996 explored the effect of 
a non uniform interstellar medium (density and tem- 
perature) in the cooling efficiencies within the giant 
cavities. In high density regions the cooling would be 
more effective and would produce a series of metal- 
rich cloudlets. Such cloudlets would fall back to the 
disk and produce local enhancements of metallicity. 
However, Recchi et al. (2001; 2008) presented argu- 
ments against this scenario, pointing out that ther- 
mal conduction and the hydrodynamical drag are 
likely to disrupt the condensations before they reach 
the disk, and thus the metals would be spread more 
or less uniformly throughout the galaxy. 

In the present paper we study the effects of young 
stellar clusters in the nuclear starburst (SB) region 
on the ISM of dwarf disk-galaxies. In particular, 
we calculate the amount of enhanced metallicity ma- 
terial that will end up in the intergalactic medium 
(IGM) and/or in the outer regions of the galactic 
disk. The paper is organized as follows: in §2 we de- 
scribe the ingredients of the galactic models that we 
use, in §3 we discuss the possible fate of the mass and 
metals in such galaxies. In §4 we describe a series of 
3D simulations of the interaction of SSCs with the 
galactic ISM. The results are presented in §5, and 
our conclusions in §6. 



2. GALACTIC MODEL INGREDIENTS 

Each galaxy in our study consists of a dark mat- 
ter halo, and a rotationally supported disk of gas 
and stars. The models are constructed following 
the approach described in Springel et al. (2005, see 
also Hernquist 1993 and Springel 2000). A near- 
equilibrium galaxy model is constructed, with a dark 
matter halo that follows a Hernquist (1990) pro- 
file, and a disk with exponential surface gas density. 
Then, using an iterative procedure, the vertical gas 
profile is determined self-consistently for a particular 
effective equation of state. 
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2.1. The halos 

We modeled the dark matter mass distribution 
with the Hernquist profile (Hernquist 1990), 
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where r is the radial distance, and Mdm is the dark 
matter mass. Two remarks can be made about 
the shape of this profile: i) in the inner parts of 
the galaxy it agrees with the profile given by the 
Navarro, Frenk & White (1996, henceforth NEW) 
fitting formula, ii) Its faster than decline makes 
the total mass to converge, allowing the construc- 
tion of isolated halos without the need of an ad-hoc 
truncation. The a parameter is related to the scale 
radius of the NEW profile by: 



a ^ r,V2[Zn(l + c) -c/(l + c)] 



(2) 



where c = r2oo/?'s is the concentration parameter, 
r2oo is the radius at which the enclosed dark matter 
mean density is 200 times the critical value (where 
the critical density is pcrit = 3iJ^/87rG). We have 
set c = 9. 

2.2. The disk 

We modeled gas and star disk components, with 
an exponential surface density profile (in the radial 
direction of the galactic disk, R) the scale of length 

^aiR) = T% exp{-R/Ro), (3) 



27ri?2 
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cxp(-i?/i?o), 



(4) 



so that the total mass of the disk is Md = {Mg + 
M^) = mdMtot, where is a dimensionless param- 
eter (fixed in this work to nid = 0.041), and Mfot 
is the total mass of the galaxy, incuding the dark 
matter halo (i.e. Mtot = Md + Mh, where Mh is the 
mass of the halo) . The vertical mass distribution of 
the stars in the disk is specified by the profile of an 
isothermal slab with a constant scale height H. The 
3D stellar density in the disk is hence given by. 



p*(i?, z) 



_^sech^ (_f_) exp 



R 
Ro 



(5) 



In all models, the H/Rq was fixed to 0.2 whereas the 
relative content of gas in the disk, the rest of the disk 
mass is in stars, was set to 0.35. On the other hand, 
the disk spin fraction (fraction of the total angular 
momentum in the disk), was set, as it is customary, 
to the value of m^. These values are typical of low- 
mass galaxies (e.g., Springel et al. 2005). 



2.3. Vertical structure of the gas disk 

The vertical structure of the disk at the begin- 
ning of the simulation is computed self-consistcntly 
with the halo and stellar disk components as follows. 
At t = the gaseous disk has uniform temperature 
(fixed to 1000 K). Consider the vertical structure of 
the gas disk is in hydrostatic equilibrium: 



dz 



"fP dz 



(6) 



where, $t = ^h+^d is the total gravitational poten- 
tial ($/i and ^d are the halo and disk gravitational 
potentials, respectively). For a given $t, the solu- 
tion of this equation is constrained by the condition 



^s(^) = / Pg{R-,z)dz, 



(7) 



where Eg (i?) is the gaseous surface mass density (see 
also eq. [3]). One can obtain the vertical distribution 
by integrating eq. ([6]), with a central density value 
at the midplane (z = 0) and an effective equation of 
state of the form P ~ P{p), see Springel et al. 2005.) 
The density starting value is guessed, then adjusted 
in an iterative scheme until the desired surface den- 
sity (eq. [3]) is recovered. This process is repeated for 
different radii to obtain an axisymmetric gas density 
distribution. 

2.3.1. Cooling 

For the non-adiabatic cases, we adopt the cooling 
functions computed in coUisional ionization equilib- 
rium (CIE) by Sutherland & Dopita (1993). With 
the energy loss per unit mass of the gas by. 



du\ 

Tt)^ 



^netiP,T) 



(8) 



We must note that the cooling function we used 
in our models assumes solar abundances. However, 
the kind of galaxies treated in this paper should be 
less metallic and thus the cooling should be lower. 
For the less metallic galaxies the results would shift 
from our radiative models to the adiabatic ones. 

3. THE FATE OF THE GAS INJECTED BY 
YOUNG STELLAR CLUSTERS 

Depending on the total amount of kinetic energy 
injected by the SSCs two outcomes are possible: if 
the kinetic energy is larger than the gravitational 
potential of the galaxy, all material will escape and 
integrate to the IGM; otherwise, the gas will return 
to the disk. If the material returns and reincorpo- 
rates to the galaxy two things can happen: it can 



4 



RODRIGUEZ-GONZALEZ, A. ET AL. 



fall back close to the galactic center or far from it. 
Considering that GWs are confined into a cone, ma- 
terial that does not reach a high altitude will fall 
close to the galactic center, while gas that reaches 
a high altitude can fall back at large galactoccntric 
radii. 



3.1. Mass loss 



measures 



In order to quantify the mass loss in numerical 
models one has to determine which material will be 
considered lost to the IGM. A natural estimation 
can be derived by comparing the velocity of a parcel 
of gas and the escape velocity at its current posi- 
tion. This has been done for instance in the previ- 
ous works of MacLow & Fcrrara (1999) and D'Ercole 
& Brighenti (1999). This works, however, used Eu- 
lerian codes in which material that leaves the com- 
putational domain can not be further analyzed, and 
therefore was automatically considered unbound. In 
reality, some of this material would return to the 
galaxy. Our models, that use the SPH (Smooth 
Particle Hydrodynamics) integration scheme, are not 
limited by the extent of the computational domain, 
they thus were able to trace all the material for the 
duration of the simulation. We will consider the fol- 
lowing more stringent criterion: only material that 
has left a cylindrical region of R„iax~'20 kpc and 
^maa:=|20| kpc (scc also D'Ercolc & Brighenti 1999) 
and has a velocity larger than the escape velocity 
at its current position will be considered effectively 
lost. 

Mass loss in our models are thus carried out by 
first computing the total mass of unbound (ejected 
from the galaxy) gas M^j. To do this, we consider 
the total galaxy gravitational potential given by, 



$-r(r) = $,,(r) -f $d(r), 



(9) 



where r is the distance to the center, and $ft and 
$d are the dark matter halo and disk components, 
respectively. Gas particles outside the cilinder with 
velocity larger than the escape velocity 



(10) 



are considered to be unbound. We define a 'total 
mass ejection efficiency' as the ratio of the mass 
ejected (unbound) to the total gas mass. 



Mn 



(11) 



The mass ejection efficiency has two limiting cases: 
6n when all the galactic ISM is bound, and 



Cm — > 1 when all the galactic gas is expelled from 
the galaxy. 

One can also define a "metal ejection efficiency" 

as 



Mr 



(12) 



where Mc is the total mass injected by the SSCs, 
and Mc,ej the mass, also originated in the SSCs that 
is unbound. The mass injected by a star cluster 
contains the metals produced in the lifetime of the 
massive stars, t ^ AQ Myr. For that time-span, us- 
ing the metal yields calculated by Meynct & Maeder 
(2002), Tenorio-Tagle et al. (2005) showed that the 
metallicity of the combined stellar winds remains at 
a value of ^ 0.5 Zq from the start of the SB, up to a 
t ~ Z Myr evolutionary time, and then rapidly grow 
to ^ 15 Zq at i ss 7 Myr, to decrease gradually to 
- 3 at t w 20 Myr. 

The metal ejection efficiency has two limiting 
cases: — >■ 0, when all the metallic gas injected 
by the nuclear starburst is kept in the galaxy disk, 
and ^2 — > 1 case, when all the metallic gas processed 
in the burst is expelled from the galaxy. Galaxies 
with ^2 1 would be metal poor, while the IGM 
around them should be metal rich. On the other 
hand, metallic gas injected from the SBs region with 
Vg < Vesc is going to return to the galactic plane, 
and the ISM of such galaxies will be contaminated 
with the metals produced in the nuclear starburst. 

3.2. Disk contamination measures 

As mentioned above, bound gas ejected by the 
SSCs will return to the galactic disk, either close 
to the center if it never reached a high altitude, or 
spread to large galactoccntric radii if it did. This 
mechanism is of course a continuum, and "high" or 
"low" altitude are just relative terms. To quantify 
the contamination of the disk by material originated 
in the SSCs we will consider as 'high altitude' a dis- 
tance of at least 3 times the scale height of the disk, 
H. 

One can define a new parameter, the fraction of 
bound gas (that is inside a cylinder of radius Umax 
and with velocity less than the escape velocity) that 
reaches at least an altitude of 3H: 



Xb 



(13) 



where M^j^h is the total bound mass (including the 
ISM and the material from the SSCs) that reaches 
an altitude of > 3H. This mass will be rein- 
serted to the galactic disk at different radii. How- 
ever, since it includes material both from the SSCs 
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(with high metalhcity) and from the ISM around the 
SSCs (lower metallicity) it is not clear how much 
the metallicity of the disk will be enhanced by this 
mixture. For that reason we define the fraction of 
metallic bound gas that leaves the galactic disk and 
reaches an altitude between 3i/ and 'Lmax as, 



Xb,: 



c,3H 



Mr 



(14) 



where Mc^sh is the bound mass injected by the stel- 
lar cluster that reaches an altitude 3H < z < Zmax, 
and Mc the total mass injected by the SSC wind. 
This gives a better measure of the amount of high 
metallicity material available to spread throughout 
the disk. Nevertheless, a detailed study of how it will 
be distributed in the disk is necessary for a complete 
disk contamination model, which is beyond the scope 
of this paper. 

4. NUMERICAL SIMULATIONS 
We have performed adiabatic and radiative 3D 
N-Body/SPH simulations that model the effect of 
a central compact starburst, in which mechanical 
energy is injected by stellar cluster winds, using 
GADGET-23. All simulations were done with 3 
types of particles: disk, halo, and gas particles, and, 
unless otherwise stated, each consist of 30000, 40000, 
and 30000 particles, respectively. Two separated 
simulations were run for each model, one adiabatic 
and a second one including radiative losses. 

Most simulations (low-resolution simulations) 
use a timestep ^ 1 — 3 x 10^ yr and a Plummer 
softening of 0.25 kpc, 0.5 kpc, and 0.25 kpc, for gas, 
halo, and disk particles, respectively. We ran a grid 
of models using different galaxies masses and star- 
burst energies. The parameters of the models are 
described below. 

4.1. The host galaxies 

We constructed ten isolated galaxies as described 
in Springel (2000) and in §2 of this paper. The galax- 
ies are named from Gl to GIO, they have ISM masses 
in the range of 5 X 10^ - 10^° Mq, similar to those 
used by MacLow & Ferrara (1999), the total masses 
are in the 4 x 10^ - 1 x 10^^ Mq range. The disk 
masses of the galaxies are given by Md = 0.041A'/tot, 
and the total mass of the gas is of 35 % of the galactic 
visible mass in all of the models. 

Table 1 shows the mass of the gas and the halo, 
the circular velocity at r2oo and Rq for all the galactic 
models. 



3GADGET-2 is a code for cosmological N-Body/SPH- 
simulation on serial workstations or massively parallel com- 



TABLE 1 
THE GALACTIC MODELS. 
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4.2. The starburst regions 

In order to reproduce the observed properties of 
starburst galaxies with young stellar clusters, includ- 
ing M82, NGC253 and NGC 1569, we use a wide 
range of SBs masses. We should note that the num- 
ber of these energetic stellar asociations in the host 
galaxies were detected at optical wavelengths (e.g. 
for M82 see Melo et al. 2005; for NGC 1569 see 
Westmoquette et al. 2008). Since there is high a 
level of obscuration in the nuclear regions of M82, 
NGC 1569, and NGC 253, the real number of clusters 
would be larger. Melo et al. (2005) and Smith et al. 
(2007) have shown that in the nuclear SBs region of 
M82 the SSCs are confined to the central ~ 200 pc 
of the galaxies. 

We estimated the mass of the starburst using 
Starburst 99 models (Leitherer & Heckman 1995 and 
Leitherer et al. 199^ with the appropriate me- 
chanical energy and a Salpeter initial mass function 
(IMF) that includes stars between 1 and 100 Mq. 
In order to obtain a better estimation of the total 
starburst masses, we have calculated the mass of the 
stars between 0.1 and 1 Mq using the Kroupa et al. 
(1993) IMF. We assumed that all of the energy pro- 
duced in the SB lifetime is injected instantaneously, 
at a given time t^, inside a spherical region of radius 
Rc at the center of the disk galaxy as thermal energy. 
At this time (corresponding to ~ 10 orbital periods 
of the galactic disk at a radius of 150 pc), the gas 
in the central regions of the galaxy has already de- 
veloped a clumpy structure. Actually, in the present 
version of the code wc cannot consider the case of an 
injection that extends over a finite time period. 



putcrs with distributed memory (Springel et al. 
Springel 2005) 



2001 & 



*http: / / www. stsci . edu / science / starburst99 / 
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Fig. 1. Time sequence of cuts of number density in xz-midplane (y = 0) for the model S8+G6 (see Tables 1&2). The 
integration times in years are given in the upper left of the frames. The top panels (a, b and c) are the adiabatic models, 
and the bottom panels are the radiative ones. 
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TABLE 2 
THE STARBURST MODELS. 
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The total injected energy is obtained from 
ftf 

E,^ / Lra{t)dt, (15) 

where t/ is the typical SBs lifetime, and Lm(t) the 
mechanical luminosity as a function of time. The 
most important contribution to the mechanical lu- 
minosity happens at ^ 40 Myr (see also Leitherer et 
al. 1999) and after that time the mechanical lumi- 
nosity decreases drastically, due to the fact that the 
most massive stars have died by then. Table 2 shows 
the total energy and the respective mass of all the 
starburst models (Sl-Sll). 

4.3. The grid of models 

We computed a total of 220 numerical simula- 
tions of GWs produced by the nuclear starburst 
events, with the different total masses of the galaxies 
Gl-GlO (see Table 1), combined with the different 
starburst energies Sl-Sll (see Table 2), all adiabatic 
as well radiative (110 models each). The total en- 
ergy injected by the starburst of the models Sl-Sll 
was distributed in the gas particles inside a sphere 
of Rc = 150 pc at the center of each of the galactic 
models Gl-GlO. 

We choose this size for the region in which the 
energy of the starburst is injected, because it is rep- 
resentative of the sizes of observed starbursts. For 
example, in M 82 the starburst region has a i?c ~ 175 
pc radius (Moorwood 1996), and in NGC 1569 A the 
starburst hole has & Rc ^ 100 pc radius (Fujita et al. 
2009). An exploration of the effects on our results of 
varying Rc is shown in Section 6. 

The number of energetic particles {N^^ and iVJp, 
for the adiabatic and radiative cases, respectively) 




7.0 7.5 8.0 8.5 9.0 9.5 10.0 
log,„ (M, [Mj) 

Fig. 2. Mass {£,m, top) and metal ejection efficiency (5^, 
bottom) countours as a function of total mass of the gas 
and total starburst masses, for adiabatic models (upper 
and lower panel, respectively). 

is related to the total number of particles of gas in 
the initial conditions. The energy per particle in 
the starburst region is Ei = Ec/Ncp, the number of 
energetic points used in the galactic models, is shown 
in Table 1. All the cluster particles (the energetic 
particles) are tagged in order to estimate the fate of 
the metals produced during the starburst lifetime. 

5. RESULTS 

Wc have let the galaxy models evolve on their 
own for tc = 1.5 x 10'' yr. At the time tc, the to- 
tal energy of the SB is injected in the center of the 
galaxy. MacLow & Ferrara (1999) injected the en- 
ergy from the cluster in a continuous way over a pe- 
riod of 40 Myr. Given the fact that the flow induced 
by the cluster takes ~ 0.5 Gyr to travel a distance 
of ~ 10 kpc (i. e., the size of the galaxy), the dif- 
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7.0 7.5 8.0 8.5 9.0 9.5 10.0 
log,„ (M, [Mj) 

Fig. 3. Same as Figure 2, but for the radiative models. 

ferences that might result from releasing the cluster 
energy over ^ 40 Myr=0.04 Gyr (rather than instan- 
taneously) cannot be very important. 

The ejection efficiencies were computed at t = 
1 Gyr in all our models. The values of the efficien- 
cies reach a constant value of integration times of 
~0.2-0.5 Gyr and ^0.4-0.7 Gyr for the adiabatic and 
radiative models (respectively). 

5.1. Mass and metal ejection efficiencies 

Total mass and metal ejection efficiencies have 
been calculated using the definitions given in sec- 
tion |3l In particular, the kinetic and gravitational 
potential energies as a function of the radial dis- 
tance ri (with respect to the center of the potential) 
were computed for each gas particle. A particle i 
would escape the galactic gravitational potential if 
Vi > Vesc = 2$i, where Vi is the velocity and <l?i 
the total potential at the position of the particle. 



We also traced the high metallicity gas particles 
(injected by the SBs) and low metallicity gas (galac- 
tic ISM gas). The high metalhcity gas was tagged 
with a negative passive scalar, the ISM with a posi- 
tive one. Therefore, was calculated using the gas 
particles with a negative passive scalar, while was 
calculated with all the gas particles (regardless of the 
sign of the passive scalar). 

Figure [1] shows the temporal evolution of the 
ejected gas from a central starburst with properties 
that correspond to starburst model S8 at the cen- 
ter of galaxy model G6 (see Tables 1 & 2). The 
three timesteps shown are 2 x 10^, 6 x 10^ and 10^ 
yr, for the adiabatic case -panels (a), (b) and (c), 
respectively-, and the radiative one ~(d), (e) and 
(f), respectively-. We can see the ejection of gas 
from the plane of the dwarf galaxy with a moder- 
ate ISM mass (3 x 10^ Mq), produced by a central 
starburst that ejects a mass of 1.5 x 10^ Mq. Such 
a SB corresponds to a region with ^ 9000 massive 
stars, injecting a total energy during its lifetime of 
around 5 x 10^^ erg (see the starburst models of Lci- 
therer & Heckman 1995 and Leitherer et al. 1999). 
From the Figure we see right away that the amount 
of gravitational unbound gas of the galaxy is larger 
in the adiabatic than in the radiative one. However, 
in both cases an important part of the disk gas is 
pushed above the galactic plane, and in the radia- 
tive case a significant part of this gas is returned to 
the galactic disk. 

For this model (S8-I-G6) the totality of the galaxy 
gas is pushed above 3H, however, only 43% acquires 
enough energy to leave the galactic potential and 
integrate with the IGM, while the remaining 57% 
would eventually return to the disk of the galaxy. 
However, this ISM material is gravitationally bound 
and is going to return to the galactic plane. For 
the radiative case, the percentage of the gas that is 
gravitationally unbound is of 14%, while 4.7% of the 
ejected gas gets above 3H and returns to the galactic 
disk. On the other hand, in the adiabatic case 97% 
of the metallic gas ejected by the nuclear starburst 
escapes from the galactic gravitational potential. In 
other words, a galaxy with those parameters would 
be contaminated only with 5% of the metals pro- 
duced in the massive SB. For the radiative S8-t-G6 
model, 79% of the metals escape from the galaxy. 

The mass ejection efficiencies arc presented in 
Figures [5] and [H for the adiabatic and the radia- 
tive model grids, respectively. In these figures, the 
abscissa sweeps the galactic models (Gl-GlO), and 
the ordinate sweeps the SB models (Sl-Sll). 

Figure [2] shows isocontour maps of the mass ejec- 
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tion efficiency (top panel) and the metal ejection 
efficiency (bottom panel) for the adiabatic mod- 
els as a function of the SB mass and the total gas 
mass of the galaxy. For the mass ejection efficiency 
(top panel) we plotted isocontours for = 0.01, 
0.10, 0.3, 0.5, 0.7, 0.9 and 0.99. In this figure, we 
can see that for the more massive galaxies (in which 
the gas mass is larger than lO^M©) the fraction of 
unbound gas is less than ^ 0.25, even for a rather 
large nuclear starburst. However, galaxies with an 
ISM of mass lower than lO^M© will lose virtually all 
of their gas if they undergo a nuclear starburst that 
injects more than IO^Mq. 

From the = 0.01, 0.10, 0.3, 0.5, 0.7, 0.8, 0.9 
and 0.99 isocontours in the bottom panel one can see 
that most of the gas ejected from the cluster escapes 
the galaxy if the starburst mass is > IO^'^Mq. Also, 
most of the ejected gas escapes if the galaxy has a 
total gas mass < lO^M©. On the other hand, for 
starburst masses < lO^'^M© and galactic gas masses 
> 1O*M0 most of the ejected metals are retained by 
the galaxy. From the bottom panel of Figure [2] we 
also see that a starburst with a mass > 10^'^Mq will 
eject more than 50% of its wind into the IGM. 

Figure |3] shows the (mass ejection) and 
(metal ejection) efficiencies as a function of galactic 
gas mass and starburst mass computed from the ra- 
diative models. Both the and contour maps 
(top and bottom panels of figure [3]) are qualita- 
tively similar to the ones computed from the adia- 
batic models (see figure[2]) . Even though the detailed 
shapes of the contours vary, the main difference be- 
tween the radiative and adiabatic cases is the follow- 
ing: in order to obtain similar ejection efficiencies, 
one needs a starburst mass larger by a factor of ^ 2 
in the radiative case. 

5.2. Material returned to the galactic disk 

Figure |4] shows the efficiencies of gas {xb, top 
panel) and metal mass (xb,z, bottom panel) that 
reaches at least 3 scale heights (H) above the disk 
but remains gravitationally bound, computed from 
the adiabatic models. Both panels show two low 
efficiency regions divided by a high efficiency ridge 
(which runs diagonally from low to high SB and 
galactic gas masses). In the high galactic mass/low 
SB mass region, the SB ejection mostly remains close 
to the galactic plane, and docs not reach the re- 
quired 3H height above the disk plane (in order to 
be counted in the Xb and Xb,z efficiencies). In the low 
galactic mass/high SB mass region, the SB ejection 
mostly escapes from the galaxy. From the Xb,z dis- 
tribution (bottom panel), we see that models with 




7.0 7.5 8.0 8.5 9.0 9.5 10.0 
log,„ (M, [Mj) 

Fig. 4. Fractions of bound gas (xt, top) and bound 
metallic gas {Xb,z) that reacli an altitude of at least 3H, 
for the adiabatic models. 
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7.0 7.5 8.0 8.5 9.0 9.5 10.0 
log,„ (M, [Mj) 

Fig. 5. Same as Figure 4, but for the radiative models. 



Mg - IQS ^ lO^Mo and Msb ~ 10^ ^ IO^Mq are 
able to redistribute more than ^ 50% of the ejected 
metals over the galactic ISM. 

Figure [5] shows the Xb (top panel) and Xb.z (bot- 
tom panel) efficiencies computed from the radiative 
models. The efficiencies show the same qualitative 
behaviour (with two low efficiency regions separated 
by a high efficiency ridge) as the ones calculated from 
the adiabatic models. However, the maximum values 
of the efficiencies are lower. In the adiabatic case, the 
maximum values of the gas mass redistribution effi- 
ciency is of ~ 20% while a maximum value of only 
~ 10% is found for the radiative models. For the 
metal redistribution efficiency, values above ^ 70% 
are obtained for the adiabatic models, while maxi- 
mum values of only ^ 20% arc obtained from the 
radiative models. 

In Figure 6, we show the distribution of metals 
that have fallen back onto the galactic disk from the 
starburst ejection of model S7-I-G9, after a IGyr time 
evolution. The SPH particles corresponding to the 
collapsed material form a distribution with a height 
of ^ 300 pc (about the initial scale height of the 
galactic disk), and with a radial extent of ~ 2 kpc 
(see Figure 6), which is about twice the size of the 
initial scale of length (Rq) of the G9 galaxy (see Ta- 
ble 1). This is an illustration of the fact that in all 
of our models the infalling metals are redistributed 
over a large region of the galactic plane. 

6. THE EFFECTS OF VARYING THE RADIUS 
OF THE INJECTION REGION AND THE 
RESOLUTION OF THE INITIAL 
STARBURST 
In all of the models presented above, we have as- 
sumed that the energy from the starburst is injected 
in a spherical region of radius Rc = 150 pc (which is 
a representative value of the observed sizes of star- 
burst regions, see §4.3). In this section, we explore 
the effects on the ejection efficiencies of varying the 
value of Rc- 

We focus on three of the starburst-|-host galaxy 
combinations which we have studied above (see 
Tables 1 and 2). We choose models S5-I-G3 (a 
small galaxy with an intermediate energy starburst, 
S10-I-G7 (an intermediate mass galaxy with a strong 
starburst) and S11-I-G9 (a high mass galaxy with a 
strong starburst). For these three models, we com- 
pute radiative simulations in which the starburst en- 
ergies are injected in spheres of radii Rc ~ 50 to 
150 pc (in 20 pc steps). The number of SPH parti- 
cles within the spheres scales as the volume of the 
initial starburst (45 particles for Rc = 50 pc to 550 
particles for i?c = 150 pc). 
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Fig. 6. Distribution of the SPH particles ejected by the central starburst for model S7+G9 after IGyr. The inner region 
of the domain is shown, and the location of the SPH particles give the spatial distribution of the metals that have been 
deposited back on the galactic disk. The left plot shows the distribution of positions projected onto the galactic disk, 
and the right plot shows the projection of the positions of the SPH particles projected onto a plane which includes the 
rotation axis of the galactic disk. The greyscale shows the distribution of the mass density on these two planes. 

In Figure 7, we show the mass ejection effi- 
ciency (see Eq. [TT]) computed for models S5+G3, 
S10+G7 and S11+G9 as a function of the initial ra- 
dius i?c of the starburst. We see that there is clear 
dependence of on i?^; with differences by fac- 
tors of ^ 1.7 models Sll-hG9 and S10+G7, and over 
- 3.5 for model S5+G3. 



We should note that the models with the smaller 
energy injection regions (with Rc — 50 pc) have 
smaller mass ejection efficiencies than the models 
with larger energy injection regions (with Rc = 
70 pc, see Figure 7). This is probably due to the 
fact that in the smaller volumes (covering the denser, 
inner region of the galaxy) the radiative cooling is 
stronger. In the models with even larger energy in- 
jection regions (with Rc > 90 pc, see Figure 7), the 
momentum per unit mass given to the galactic gas 
is lower, and therefore less material succeeds to es- 
cape from the galactic potential. From these results, 
we conclude that in our intermediate and high mass 
galaxy models, variations of factors of ~ 2 in the ef- 
ficiencies are to be expected for starbursts (with the 
same energies) with initial radii in the Rc = 50 to 
150 pc range. 

We have also explored the effect of changing the 
total number of SPH particles in the simulations. We 
again choose models S5-I-G3, S10-HG7 and S11-HG9 
(see above and Tables 1 and 2) as representative of 
our model grid, keeping a constant, Rc = 150 pc 




80 100 120 

Cluster Radius (R.) 



Fig. 7. The mass ejection efficiency computed for 
models S5-fG3, S10H-G7 and S1H-G9 as a function of 
the initial radius Rc of the starburst. 
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initial radius for the starbursts, and compute radia- 
tive simulations with more than one order of mag- 
intud more of gas, halo and disk particles (1500000, 
2000000 and 1500000, respectively) than the previ- 
ously computed models (see §4). In these new runs 
we have also changed the softening length to 50 pc, 
for all the components. In all the new models the 
mass ejection efficiency increases by a factor of the 
order of 2. For model S5-I-G3 this change is a factor 
of 1.5, while S104-G7 and S11-I-G9 suffer changes 
of factors of ~ 2.21 and ^ 1.6, respectively. 

Therefore, we seem to obtain at least a partial 
convergence for the values of the ejection efficiencies 
in simulations with ^ 10^ gas particles. However, 
the computing time necessary for these simulations 
is quite large (~ 30 times longer than the 30000 gas 
particle models), so that our parameter study has 
been done at lower resolution. 

7. DISCUSSION 

MacLow & Ferrara (1999) computed radiative, 
axisymmetric simulations of the ejection of gas from 
galaxies due to central starbursts. Their models dif- 
fer from ours in that the starburst energy is injected 
as a continuous source term during the 50 Myr of 
their simulations, while in our models the energy of 
the starburst is injected instantaneously at the be- 
ginning of the simulations. 

If we convert their starburst mechanical lumi- 
nosities into total energies (by multiplying the lumi- 
nosities by 50 Myr), and using the relation between 
starburst energy and mass of Table 2, we can find 
correspondences between their models and ours. We 
then find that both the mass and metal ejection effi- 
ciencies computed by MacLow & Ferrara (1999) are 
generally consistent with the ones that we have com- 
puted, with typical differences of factors ~ 3. We 
judge this to be a satisfying agreement given the dif- 
ferences between our models and the ones of MacLow 
& Ferrara (1999). 

Interestingly, the highest SB mechanical en- 
ergy/lower galactic mass model of MacLow & 
Ferrara (1999) approximately corresponds to our 
S5-hG2 model (with M^s = 3 x 10* Mq and Mg = 
3 X lO"^ M0). As can be seen in Figure 3, this model 
lies within the transition from low to high efficiency 
regimes. Therefore, the models of MacLow & Ferrara 
(1999) mostly missed the transition between low and 
high mass loss efficiencies which we are describing in 
the present paper. 

D'Ercole & Brighcnti (1999) computed a set of 
5 models of galaxies with Mg < 1.4 x lO^M©, sim- 
ilar to our galaxy models G1-G5, in which they in- 
ject energies in the range of our S6-S8 starbursts. 



Their "PEXT" model (see Table 2 of D'Ercole & 
Brighenti 1999) has parameters close to our S8-I-G4 
model (see Tables 1 and 2). The PEXT model has 
a Cm = 0.31 mass ejection efficiency, which is consis- 
tent with the Cm = 0.45 efficiency obtained from our 
S8+G5 model. D'Ercole & Brighenti (1999) have 
studied different possibilities (c. g., having galaxies 
with no rotation, or having fixed dark matter mass 
distributions) which have not been explored in our 
present work. 

Fragile et al. (2004) computed mass and metal 
ejection efficiencies produced in galaxies with IO^Mq 
(similar to our galaxy model G9) in which they inyect 
energies in the range of our model S3-S7 starbursts. 
Their model 5 (see Table 1 of Fragile et al. 2004) has 
parameters close to our S7+G9 model. Model 5 has 
a Cm = 0.016 mass ejection efficiency and Cz = 0.99 
metal ejection efficiency (see Table 2 of Fragile et al. 
2004), which is consistent with the Cm = 0.006 and 
Cz = 0.805 efficiencies of our model S7-I-G9. Frag- 
ile et al. (2004) have explored the metal ejection 
efficiencies for supernova events at different radial 
positions in the galactic disk. They find that a sub- 
stantial fraction of the metals are retained in their off 
center starburst models. These result shows that if 
we consider off-galactic center starbursts, the mass 
and metal ejection efficiencies will be substantially 
lower than the ones that we have computed in our 
present paper. 

It is worth to mention that the strong winds pro- 
duced by the combination of individual stellar winds 
and SNc explosions not only produce the GWs, but 
also have important consequences on the star for- 
mation history of their host galaxy as well. They 
are the most efficient way to redistribute metals to 
the ISM (and/or the IGM). This produces a feed- 
back mechanism that facilitates the gravitational col- 
lapse of the ISM to form stars. The material with 
higher metallicity cools more efficiently, and it there- 
fore forms high density condensations on shorter 
timescales than metal poor gas (Scannapiecco et al. 
2005, Sutherland & Dopita 1993). At the same time, 
the energy feedback from these strong winds can heat 
and disrupt cold gas clouds, thus inhibiting star for- 
mation. 

Observations and theory suggest (e.g. Larson 
1974, White & Rees 1978, White & Frenk 1991) that 
large outflows should be able to develop in small sys- 
tems because of their shallower potential wells. This 
is because, according to hierarchical galaxy forma- 
tion scenarios, large systems are formed by the ag- 
gregation of smaller ones. Thus, the energy feedback 
from strong winds should have a more important ef- 
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feet for systems in the early stages of formation. 

Our models provide quantitative predictions of 
the amount of metals which are ejected into the IGM 
and which are redistributed within the host galaxy. 
These predictions could be used in a direct way for 
computing models of the star formation history of 
the host galaxy, and could also be included in models 
of the history of galaxy formation. 

8. CONCLUSIONS 

We computed 220 numerical simulations of GWs 
produced by nuclear starburst energy injected inside 
a spherical region (of radius Rc) at the center of 
the disk galaxy. We used different total masses of 
the galaxies Gl-GlO, combined with different star- 
burst energies Sl-Sll, all adiabatic as well as radia- 
tive (110 models each). We have performed adia- 
batic and radiative 3D N-Body/SPH simulations of 
GW models using GADGET-2. the radiative models 
have been calculated using a solar metallicity cooling 
function. Finally, the adiabatic models were calcu- 
lated in order to obtain the maximum possible value 
of the mass and metal ejection cffciency in our mod- 
els of galactic winds, and do not necessarily corre- 
spond to real cases of starbursts in galaxies. 

From the grid of models we obtained the mass 
and metal ejection efficiencies as a function of the 
mass of the gas in the galaxy and the total mass of 
the nuclear starburst. 

In the same way, we computed the efficiencies of 
gas and metal mass that reaches at least 3 H above 
the disk but remains gravitationally bound for both 
adiabatic and radiative regimes. We present these 
efficiencies as a function of the mass of interstellar 
medium in the galaxies versus total mass of the nu- 
clear starburst. 

We show that for compact starburst regions (e.g. 
inside the SSC) with masses larger than 10^ Mq in 
the adiabatic case, the metal injected by the nuclear 
starburst is expelled to the IGM (at least 90% of 
this gas). We show that in galaxies with gas masses 
between 10^ to 10^° Mg and with a massive nuclear 
starbursts the available metals are freed from the 
gravitational potential of the galaxy. In the radia- 
tive case the critical SB mass at which total loss of 
metals occurs is around 2 or 3 times larger than the 
adiabiatic one. 

In addition, we have calculated the efiiciency of 
the metal redistribution in the disk of galaxies. We 
find that the metal redistribution is more efiicient 
for intermediate and low mass starburst masses (< 
lO'^-^Mg) and galactic gas masses > IO^Mq (for both 
the adiabatic and radiative cases). 



In this paper we presented models in which we 
do not consider many elements that are likely to be 
important in starburst galaxies, among these: 

• dumpiness in the medium, 

• the generation of new stars, 

• a different fraction of gas in the galaxy disk, 

• the effect of galactic mass loading (see Fujita et 
al. 2009). 

Therefore, there is a large range of possibilities 
for future studies including differential and/or selec- 
tive galactic winds (also, Pilyugin 1993 and Recchi 
ct al. 2001, 2002). 

Finally, in our paper we have presented a limited 
resolution study (in section 6). Our evaluation of 
the dependence of the computed efficiencies on reso- 
lution is limited to a range of ~3 in resolution. This 
range is of course not very impressive, and questions 
remain as to the results that would be obtained at 
considerably higher resolutions (presently not possi- 
ble, at least for calculating relatively extended grids 
of models) . When such higher resolution models do 
become available, one will probably find different nu- 
merical values for the efficiencies, but one would ex- 
pect that the transition between low and high ejec- 
tion efficiency systems (found in the present work) 
will be preserved. 
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